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Abstract: We consider the critical behavior of two-dimensional layered Ising models 
where the exchange couplings between neighboring layers follow hierarchical sequences. 
The perturbation caused by the non-periodicity could be irrelevant, relevant or marginal. 
For marginal sequences we have performed a detailed study, which involved analytical 
and numerical calculations of different surface and bulk critical quantities in the two- 
dimensional classical as well as in the one-dimensional quantum version of the model. The 
critical exponents are found to vary continuously with the strength of the modulation, 
while close to the critical point the system is essentially anisotropic: the correlation length 
is diverging with different exponents along and perpendicular to the layers. 



PACS-numbers: 05.50.+q, 64.60.Cn, 64.60.Fr 



1 



I. Introduction 



Since the discovery of quasicrystals[l] there is a growing interest to understand their 
structure and physical properties. Theoretically it is a challenging problem to understand 
the properties of phase transitions in these quasiperiodic or more generally non-periodic 
structures. Since a non-periodic system shares some aspects with a system with quenched 
disorder, from these studies one hopes to get a better understanding about the critical 
behavior of random systems, too. 

Early studies on this field were numerical investigations on specific problems (Ising 
model[2-4], per eolation [5, 6], random walks[7] etc.) on two- and three-dimensional quasi- 
periodic lattices, whereas analytical results were obtained on layered two-dimensional 
Ising models with one-dimensional aperiodicity[8,9] or on the corresponding Ising quan- 
tum chain[10-14]. The nature of phase transitions is better understood, since Luck[15] 
generalized the Harris criterion[16] for non-periodic systems. In this relevance-irrelevance 
criterion the strength of fluctuations in the couplings in a scale of the correlation length is 
of primary importance. The perturbation caused by the non-periodicity is irrelevant (rele- 
vant) if the local energy fluctuations are smaller (greater) than the corresponding thermal 
energy. Theoretically most interesting is the borderline case, when the thermal and fluc- 
tuating energy contributions are in the same order of magnitude, i.e. the perturbation 
is marginal. In this case - according to exact results on two-dimensional layered Ising 
models[17-20] - the critical behavior is non-universal, the critical exponents are continuous 
functions of the stregth of the modulation. Furthermore, close to the critical point the 
system becomes essentially anisotropic [20], i.e. the correlation length is diverging with 
different exponents along and perpendicular to the layers. 

The Luck criterion is obtained in the frame of a linear stability analysis therefore 
its validity is restricted to such non-periodic systems where the perturbation in the local 
couplings is small. Thus the criterion is valid for quasiperiodic lattices and for such one- 
dimensional aperiodic sequences wich are generated by substitutional rules. On the other 
hand the criterion is no longer valid in its original form, if the modulation of the couplings 
follows some hiererchical sequence, such that certain couplings could become arbitrarily 
large or small. This type of hierarchical sequence[21] was introduced first by Huberman 
and Kerszberg[22] and generalized by others[23,24] to study anomalous diffusion and the 
properties of the spectrum of the Hamiltonian[25] in one-dimension. 



2 



As far as the critical properties of Ising quantum chains with a hierarchical structure 
in the couplings are concerned two conflicting results exist. From a study of the low-lying 
excitations of the system Ceccatto[26] has drawn the conclusion, that the perturbation is 
irrelevant, if the hierarchical parameter r is smaller than some critical value r c . Whereas for 
r > r c the perturbation is relevant, the Ising-type critical point of the homogeneous system 
is washed out by the perturbation and the critical behavior of the system is similar to that 
of the McCoy- Wu mo del [8], i.e. to a layered Ising model with quenched randomness. In 
contrary to this results Lin and Goda[27] obtained continuously varying surface magneti- 
sation exponent of the hierarchical quantum Ising model, thus according to this study the 
perturbation caused by this non-priodicity is marginal in the whole range of the parameter 
r. Similar conclusion is drawn from renormalization group and Monte Carlo simulation 
studies on the two-dimensional Ising model with layered hierarchical couplings [28]. 

Our aim in the present paper is to clarify this controversary and to present a compre- 
hensive picture about the critical behavior of hierarchical Ising models. For this purpose 
we consider general hierarchical sequences and investigate the condition for relevance (ir- 
relevance) of the perturbation. For marginal hierarchical modulation we perform a detailed 
study, which includes analytical and numerical calculations of different bulk and surface 
critical quantities both in the two-dimensional classical and in the one-dimensional quan- 
tum version of the model. Finally, the results are discussed in the frame of a general scaling 
theory, which is then compared to the analogous one of aperiodic systems. 

II. The model 

We consider the Ising model on the square lattice with different layered structures: the 
system is translationally invariant either along the columns (Fig la) or along the diagonals 
(Fig lb). In the first case the interaction along the layers Ki is constant, whereas it is 
modulated in the other direction and given as K 2 {k) in units of 1/ksT between neighboring 
layers at k and k + 1. In the extreme anisotropic limit K\ — > oo, K 2 (k) ^ the transfer 
matrix of the problem involves the Hamiltonian of a quantum Ising chain: 




(1) 



k 



where cr£, af. are Pauli matrices at site k and 
A fc = -2K 2 (k)/ln(tanhK 1 ). 
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In the second situation the K^(k) diagonal couplings are the same within one column 
and the quantities 

Y k = sinh[2K d (k)] (2) 

are assumed to vary in a hierarchical way. 

The hierarchical sequences we use in this paper were first introduced in economical 
problems[21] and later applied to study the so called hyperdiffusion process. For an integer 
number m the sequence 0,1,0,2, •• . is defined as: 

cifc = ar-^ k (3) 

where r is the ratio of the sequence, a is some reference value and the fk-s are natural 
numbers satisfying the relation: 

k = m fk (ml + n) , / = 0, 1,... , \i = 1,2, . . . ,m - 1 (4) 

As pointed out in Ref[29] it is possible to generalize the sequence by modifying eq(3) 
as ak = ar 9 ^ k \ where g(x) is some analytical function. Here we shall consider power 
functions: g(x) = x w , ui > 0, thus the original sequence in eq(3) corresponds to u> = 1. 

For the layered Ising models we assume the hierarchical variation in the couplings 
and in the parameters Yk, respectively. 

III. Surface magnetisation 



For the two-dimensional Ising model the surface magnetisation m s is the simplest 
order parameter, which can be most easily determined in the extreme anisotropic limit 
eq(l) through the formula: 



-1/2 

00 j 



TO.s 



3=1 k=i 



For the general hierarchical sequence containing the exponent u> the surface magnetisation 
is rewritten in the form: 

00 j 

m s = [S(X,r)]- 1/2 , S(X,r) = Y / ^ 2j r- 2nj , = £(/ fc ) w » "0 = (6) 

3=0 fc=i 



The critical coupling A c is such that [30] 



1 3 

lim - \ In Afc = 

j' — \ r>r\ 1 f ^ 



J k = l 



(7) 



and related to the hierarchical parameter r as 

A c = r - 5(w ' m) (8) 
where 5{u,m) = lim^oo rij/m- K As shown in the Appendix: 



which can be expressed in closed form for integer cus. For ui = 1, 2 and 3 it is given as: 

d(l,m) = , d(2,m) = - , d(3,m) = — —3— (10) 

m — 1 (m — l) 2 (m — l)" 3 

To calculate the surface magnetisation we do separately for uj = 1 and to ^ 1. In the 
first case, for u; = 1 one can verify that f mp = f p + 1 and / mp + M = , /ii = 1, 2, . . . , m — 1, 
from which the relations n mp = n p + p , rimp+fj, = n mp follow. Then deviding the sum 
for S(\,r) in eq(6) into m parts: 

00 m-l 00 

p=0 M=l P=0 

and using the relations between the n^-s one gets the following equation: 

1 \ -2m 

S(A,r) = £(A m r,r) — — (12) 

From this expression the surface magnetisation exponent (3 S defined by m s (t) ~ (t)^ s as 
t = 1 — (A c / A) 2 — > + can be evaluated as in Ref[18]. According to eq(6) close to the 
critical point S(t) ~ (t)~ 2/3s , and the critical exponent corresponding to eq(12) is given by: 

In 

Ps = \\~^ J (13) 
2mm 

In the special case m = 2 one recovers the result in Ref[27]. According to eq(13) for 
uj = 1 the critical behavior of the hierarchical Ising model is non-universal, since (3 8 is 
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a continuous function of the ratio r. This functional dependence is shown on Fig 2. for 
several values of the parameter m. 

Next we turn to discuss the situation for uj ^ 1, in which case the surface magnetisation 
is studied on large finite sequences. Let us consider the system of length j = m N — 1 = L— 1, 
i.e. after N generations, and define in analogy with eq(6): 



m — 1 J m —1 



S N (X,r) = 1 + J2 U X k 2 = E A- 2 ^- 2 ^ (14) 

j=l k=l j=0 

This quantity satisfies the relation: 

i y"2m 

S N+1 (X,r)= l _ x2 S N (X,r) , X = X~ L r~ nL (15) 

As shown in the Appendix in the large L limit the finite size corrections to the criticality 
condition in eq(8) are logarithmic: 

n L = 5(u,m)L ^— -iV w_1 + O (N w ~ 2 ) (16) 

m — 1 

thus the parameter X in eq(15) in leading order is given as X = (A c /A) L r^/^ m_1 ^ Ar 
and at the critical point one gets the relation: 

WAo) = 1 ~ r ^S N (X c ) , Z = r^-W- 1 (17) 

where r is related to A c through eq(8). 

For large N the parameter Z in eq(17) scales differently for uj = 1, oj < 1 and u > 1, 
respectively, and the corresponding functional form of Sn{X c ) is also different in the three 
cases. In the borderline case ui = 1, which is studied already before, Sn(X c ) has a power 
law dependence on the size of the system: 

S N (X C ) ~ (L) 2Xa , w=l (18) 

where according to finite size scaling[31] x s is the scaling dimension of surface spins. Using 
the value of 5(1, m) in eq(10) one gets from eqs(17) and (18) that x s corresponds to (3 S in 
eq(13), thus the scaling relation (3 8 = vx s is satisfied, since the correlation length critical 
exponent for the two-dimensional Ising model is v = 1. 
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By finite size scaling one may investigate the surface critical behavior of the system 
at the right end of the chain taking L = m N spins, t Denoting the inverse square of the 
finite size surface magnetisation as Sn{\ t) one can write the simple relation: 



S N (\r) = l + \-'r-^S N (\,r) (19) 

from which x s the scaling dimension of surface spins on the right end of the chain is given 
by: 



lnr ln 



i-K m 

l-A? 



s s lnm 2 ln m 



(20) 



thus x s {\ c ) = x^X' 1 ). 

For < uj < 1 Z in eq(17) goes to zero, and the relation in eq(17) in leading order 
reads as 

S N+1 (X C ) =m(l + uN aj - 1 \nr)S N (X c ) , < u < 1 (21) 

which is solved by the function Sn{X c ) ~ m N r N " . Now using from finite size scaling 
theory [31] that m N ~ L ~ one obtains with v = 1 for the temperature dependence 
of the surface magnetisation: 

ms {t)^ t V2 r -(\^t\/lo gm r/2 ^ Q<u<1 (22) 

Thus the surface magnetisation in this case vanishes with the same exponent f3 s = 1/2 
as the homogeneous model with r = 1, however there is a logarithmic correction. Con- 
sequently the hierarchical perturbation in the copulings for < u < 1 is marginally 
irrelevant. 

The behavior of the system is completely different for u > 1. In this case one has to 
study separately r > 1, when Z in eq(17) goes to infinity and r < 1, when Z goes to zero. 
In the first case the asymptotic relation 

S N+1 (\ c )=r 2 " N "- 1 S N (\ c ) (23) 

is satisfied by the function 

SW(A c )~r 2JVW (24) 



t If the chain consists of m N — 1 spins as before, then it is symmetric to its center and 
the surface magnetisation is the same at both ends. 
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Therefore the coupling dependence of the surface magnetisation is anomalous: 

m s {t) ^ r -(|iogt|/iogm)" ^ u>1 ^ r>1 ( 25 ) 

it vanishes faster than any power as A — > A c . 

On the other hand for r < 1 (and u; > 1) Z goes to one, and the asymptotic relation 

S N+1 (X C ) = (l + r^^" 1 ) 5at(A c ) (26) 

implies lini/v-^oo -SV(A C ) < oo, since the product riiv=iVo ^1 + r m ~ lN j is convergent. 
Consequently the surface magnetisation stays finite at the critical point and the phase 
transition at the surface is of first order. According to eq(26) the surface magnetisation 
approaches its limiting value m s (0) in an anomalous way: 

m s (t) -m a (0) ^ r llogVlogmr- 1 2a;/(m-l) ^ 

where again the correspondence L ~ t~ x was used. The surface magnetisation for r = .92 
and different values of ui > 1 are shown on Fig 3. 

Next we turn to present results about the critical behavior of the model at the (1,1) 
surface, when the non-periodicity in the couplings is given in the diagonal direction (Fig 
lb). The criticality condition now reads as [32]: 

1 j 

lim - Vlny fc = (28) 

J k=l 

where the parameter = smh[2Kd(k)] plays a similar role as Afc for the quantum Ising 
model in eq(7). To make use further this analogy we assign the hierarchical modulation 
in eq(3) to the generalized parameters. As a consequence the critical value Y c and the 
ratio of the sequence r are related as in eq(8) for the quantum model. 

To study the magnetisation on the (1,1) surface we use a numerical procedure by 
Hilhorst and van Leeuwen (HvL), which is originally developed for the triangular lattice 
and based on a repeated use of the star-triangle transformation [33]. From the square 
lattice on Fig lb one can obtain the triangular lattice by connecting vertically next-nearest 
neighbors with non-vanishing couplings. The numerical procedure also works for this 
triangular lattice, however, we shall only study the problem on the square lattice. Details 
of the method together with results on the surface critical behavior of different type of 
non-periodic triangular lattice Ising models will be presented elsewhere [34]. 
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The HvL method enables one to obtain very accurate numerical estimates on the 
magnetisation at the (1,1) surface as well as to study the decay of surface correlations. 
The obtained results on the surface critical behavior for marginally irrelevant (uj < 1) 
and marginally relevant (u> > 1) perturbations are qualitatively the same as we found for 
the corresponding quantum chain. In the truely marginal case - u> = 1 - the results even 
quantitatively agree with those obtained on the (1,0) surface. We are going to discuss 
these results in more detail. 

First investigating the surface magnetisation one can see that in the limit r — > 
m s (T) = 1 for T < T c and m s (T) = for T > T c , i.e. it stays finite at the criti- 
cal point. This anomalous behavior is due to the topology of the lattice in the diagonal 
direction. As r — > in the first layer (and also in any even layers) Y\ = Y c = 1/r, which 
together with K^iX) goes to infinity, resulting an ordered surface layer. For r > ac- 
cording to numerical results the surface magnetisation vanishes at the critical point. For 
large number of iterations / the critical surface magnetisation goes to zero as a power: 
m s (t = 0) ~ I -7 , where the exponent 7 > is related to the decay exponent r)\\ in eq(29) 
as 7 = ?7||/4. 

The surface magnetisation exponent (3 8 determined from the relation m s (t,r) ~ t^ s is 
found to be the same function of r as that at the (1,0) surface in eq(13). Similar conclusion 
holds for the magnetisation exponent on the right end of the chain. We obtained x s ( Y c ) = 
x s (Y~ 1 ) in close analogy with the results on the (1,0) surface. Numerical estimates on the 
surface magnetisation exponents are drawn on Fig 2. We note that the accuracy of the 
calculation is limited by the fact that the error in the magnetisation decreases very slowly 
with the number of iterations / as J -1 / 2 . 

We have also investigated the decay of critical surface spin correlations. The decay is 
given as a power: 

G s ( ril ,t = 0)^r~ vn (29) 

and according to the numerical estimates the critical exponent rj\\ can be described with 
high accuracy by the formula: 

v\\ = — — ' ^11 = — r~ ( 30 ) 

both for r < 1 and r > 1. In eq(30) f\\\ is the decay exponent on the right surface of the 
system. 

The relations in eq(30) are in conflict with ordinary scaling, which would imply r}\\ = 
2x s . They could be explained, however, within the frame of anisotropic scaling theory [35]. 
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Then the critical spin-spin correlation function at the left side of the system is expected 
to behave as: 

G s (r h t) = b- 2x °G 8 (r ll /b z ,b 1 /»t) (31) 

when lengths perpendicular to the layers transform as r± — > rj_/5, whereas scaling along 
the layers involves the dynamical exponent z: r\\ —>r\\/b z . As a consequence the correlation 
length close to the critical point diverges with different exponents in the two directions: 
~ t~ v and £y ~ t~ vz . According to eq(31) the decay exponent is given by v\\ = 2x s /z, 
thus the numerical results in eq(30) imply that the dynamical exponent is expressed by 
the sum of the two surface magnetisation exponents: 



In 



(32) 



mm 

In the following we check the validity of the anisotropic scaling hypothesis by calculating 
other critical quantities. 



IV. Other critical quantities 

A direct evidence for anisotropic scaling can be obtained from the behavior of the 
correlation length. At the critical point the perpendicular correlation length on a strip of 
width L is restricted to £j_ ~ L, thus £y ~ £j_ ~ L z . In the extreme anisotropic limit £y 
is given by the inverse gap of the critical Hamiltonian in eq(l), thus one expects for low 
laying states: 

Ei-Eo- L~ z (33) 

To check this relation we have numerically studied the excitation spectrum of the Hamilto- 
nian in eq(l) at the critical point. Using standard methods [36] a quadratic fermion Hamil- 
tonian is obtained via the Jordan- Wigner transformation, which is then diagonalized by 
a Bogoliubov transformation. Using free boundary conditions we calculated the energy of 
the first fermion modes on finite systems of size L = m N — 1 (m = 2, N = 2, . . . , 16; m = 3, 
N = 2, ... , 10). Their size dependence is found to be very accurately described by the re- 
lation in eq(33). The ^-exponents obtained through sequence extrapolation methods[37] 
coincide with the values in eq(32) up to 4 — 5 digits. 

Even more accurate estimates can be obtained if the energy of fermion modes are 
calculated in a second order perturbation expansion[15]. At the critical point their size 
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dependence follows from the relation: 



L L 
j=l k=l 



1 -1/2 



2 

j+k 



(34) 



Evaluating eq(34) up to L = 2 24 and L = 3 15 the accuracy of the estimates on the z- 
exponent has been increased up to 7-8 digits. 

Next we turn to study the behavior of the specific heat on finite systems. According 
to anisotropic scaling the singular part of the bulk free energy density in a finite system of 
size L behaves as[35]: 

f(t,L)=b-^f(b 1 /%L/b) (35) 

This expression is in accord with the finite size scaling behavior of the singular part of the 
critical bulk energy density e(L) ~ L _Xe , since according to our numerical estimates: 



Xp — z 



(36) 



From eq(35) the finite size dependence of the specific heat at the critical point is given by 
C(t = 0, L) ~ L~ a -I v , where the specific heat exponent 



a = 1 — z 



(37) 



is negative for r ^ 1. As shown in Table 1 this relation is also satisfied. 

Finally, we report on our results on the scaling behavior of the surface energy density 
at the critical point: e s (L) ~ L~ Xes . On the left surface we obtain: 



x es — z -\- 2x 6 



(38) 



and a similar expression is valid on the right surface. 
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V. Discussion 

In this paper the effect of a layered hierarchical structure on the critical properties of 
the two-dimensional Ising model is studied by analytical and accurate numerical methods. 
The perturbation is found to be marginal for a; = 1, which corresponds to the original 
Huberman-Kerszberg series. In this case the critical exponents are continuous functions 
of the hierarchical parameter r, for the whole range of < r < oo. These results are 
in accordance with the findings of Lin and Goda[27] on the surface magnetisation in the 
m = 2 model, however they are in conflict with the results of Ceccatto[26]. The failure of 
Ceccatto's calculation is due to the fact that his perturbational method is only valid for 
irrelevant inhomogeneities and can not be used for the hierarchical sequence, which is a 
marginal perturbation. Our conclusions are also consistent with the RG and MC simulation 
results of Stella et al[28] on a two dimensional hierarchically layered Ising model. 

The nature of the hierarchical perturbation is found to be different for ui < 1 and 
uj > 1. In the first case the perturbation is irrelevant, however there are logarithmic 
scaling corrections. In the other regime, for uj > 1, the perturbation is marginally relevant 
and the surface magnetisation behaves anomalously at the critical point. In the following 
we present a relevance-irrelevance criterion, which explains the above results. 

This criterion is actually a modification of the Harris criterion[16] for random systems, 
which is then generalized by Luck[15] to aperiodic systems. In this criterion one compares 
the energy Ef due to fluctuations in the couplings with the thermal excess energy E t . If the 
couplings Jk follow a one-dimensional aperiodicity and their average is J the fluctuating 
energy in a domain of size of the bulk correlation length L ~ £ behaves asymptotically as: 

L 

E fl (L) = J2(Jk-J)~L Q (39) 
fc=i 

where O < 1 is the wandering exponent characteristic to the sequence [38]. The fluctuating 
energy per spin then scales as e//(L) ~ L 51-1 . The excess thermal energy per spin is 
proportional to the reduced temperature e t {L) ~ t ~ L~ x l v . Comparing e/i(L) with e t (L) 
one assumes irrelevant perturbation for e//(L) <C e t(L), i.e. for O < 1 — v and relevant 
modulation in the opposite limit Q > 1 — v. Finally, the perturbation is marginal for 
O = 1 — z/, which condition is satisfied for the Ising model with O = 0. 

For the hierarchical perturbation considered in this paper the above criterion can not 
be used directly. The perturbation is locally not small and for r > 2 even the average 



12 



value of the coupling A is divergent. To remedy this problem we consider first the diagonal 
hierarchy in eq(2), where the parameters are given as Y k = Yir n . For a large n Y k ~ 
ex.p[2Kd(k)] ~ r n , thus the coupling itself is given as Kd(k) ~ log Yk ~ nKd(l). Therefore 
the fluctuating energy in this case is expressed in terms of the logarithms of the parameters: 

L 

^(L)~^log(y fc /F)~L n (40) 
k=i 

Similar conclusion is valid for the extreme anisotropic system, where K2(k) — > 0, K\ — > oo 
and Afc = K2(k) exp(2Ki) = Air n . To that quantum system one can assign another two- 
dimensional classical system, if K2 is kept constant and the vertical couplings K\(k) vary 
with the column index k such that their ratio remains A&. Then Ki(k) ~ log Afc ~ nKi(l), 
and the fluctuating energy is again given by: 

L 

E fl (L)^J2^g(Xk/X)^L n (41) 
k=i 

Since eqs(40) and (41) remain valid for small aperiodic perturbations they can be consid- 
ered as the general definition of the fluctuating energy both for bounded and unbounded 
sequences. 

For generalized hierarchical sequences the fluctuating energy can be obtained by 
noticing that at the critical point Ylk=i ^ k = wm ch is defined in eq(17). Since 
A = liniL^oofn^i A fc ] 1/L = 1 we obtain 

Efi(L) ~ logr 

Thus the wandering exponent O = 0, for any value of uj and a linear stability analysis 
can not decide about the relevance- irrelevance of the perturbation. In second order of the 
analysis the perturbation in the uj > 1 region (where Efi(L) logarithmically divergent) 
is predicted to be relevant, whereas for uj < 1 it is marginally irrelevant. Finally, for 
uj = 1 Efi(L) is independent of the size of the system and the perturbation is truely 
marginal. These predictions of the modified Harris criterion are in complete agreement 
with our analytical results. 

The critical behavior found at uj = 1 is very similar to that obtained in the two 
dimensional Ising model with marginally aperiodic layered interactions. In both type of 
problems the system is essentially anisotropic at the critical point and the corresponding 



UJ 



m 



1 \logm 



UJ-l 



(42) 
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z exponent is related to the surface scaling dimensions, as in eq(38). The surface spin 
correlations follow the scaling law in eq(31), whereas the thermodynamic singularities are 
consistent with the anisotropic scaling form of the free-energy density in eq(35). Finally, 
we note that some exponents, including z, x e and be obtained exactly by an 

analytical method, which can be used both for aperiodic and hierarchical marginal Ising 
chains. Details of the calculations will be presented elsewhere [39]. 

Acknowledgement: This work has been supported by the Hungarian National Re- 
search Fund under grant No OTKA TO12830. F.I. thanks for interesting discussions with 
L. Turban, B. Berche and A. Maritan. 
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Appendix 



We calculate the rij parameter in eq(6) for sequences of lengths j = m, m 2 , . . . , m N , 



n m = 1 



n m 2 = n m m + 2^-1^ 



n m N = n m N-im + iV^ — (N — 1) L 
which relations are satisfied with 



n mN =m N (l--) ^^7 + ^ 



m I * — ' m? 



From this expression the value of 5(u,m) as given in eq(9) follows. For a positive integer 
u> S(u>, m) can be expressed as: 



in) ----- <j (1 - x) [ x — 



The finite size corrections to n m N are given by: 



1 



1-x 



x=l/m 



n mN - *(u/, m)m N = - ( 1 - - ) £ J — + AT 

j=N 



= N u 



» m ) m 3 

7 j=o 



Expanding (1 + j /N) w in Taylor series and performing the summation the leading term is 
given in eq(16). 
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Figure captions 



Fig.l: Two-dimensional Ising model with layered hierarchical couplings: (a) along the col- 
umns and (b) along the diagonals. 

Fig. 2: Surface magnetisation exponent for marginally hierarchical Ising chains at the left 
surface (f3 s ) and at the right boundary (0 S ) as a function of the parameter r. The 
lines correspond to the analytical results in eqs(13) and (20) for the (1,0) surface and 
the squares are numerical estimates for the (1,1) surface. 

Fig. 3: The surface magnetisation for relevant perturbations (r = .92) at different values of 

UJ. 
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Table caption 



Table 1: Finite size estimates for the specific heat exponent of the m = 2 marginally hierarchical 
Ising model for different values of the critical coupling A c . The values in the table 
were calculated from finite size results on chains of length N, N/4 and N/16. In the 
last row the conjectured results from anisotropic scaling eq(37) are given. 
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N 



\l = 2 



XI = 3 



Xl = 4 



XI = 5 



32 

64 

128 

256 

512 

1024 

2048 

4096 



-.11269 
-.08671 
-.08008 
-.08007 
-.08154 
-.08288 
-.08378 
-.08434 



-.27620 
-.23422 
-.21749 
-.21118 
-.20886 
-.20801 
-.20771 
-.20765 



-.43191 
-.37389 
-.34636 
-.33350 
-.32748 
-.32460 
-.32319 
-.32285 



-.58306 
-.50684 
-.46609 
-.44517 
-.43454 
-.42922 
-.42679 
-.42507 



1-z 



.08496 



.20752 



.32193 



.42400 



Table 1 
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